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In normal degenerate quark matter, the exchange of dynamically screened 
| transverse gluons introduces infrared divergences in the quark self-energies 

that lead to the breakdown of the Fermi liquid description. If the core of neu- 
tron stars are composed of quark matter with a normal component, cooling by 



direct quark Urea processes may be modified by non-Fermi liquid corrections. 



We find that while the quasiparticle density of states is finite and non-zero at 
the Fermi surface, its frequency derivative diverges and results in non-Fermi 
liquid corrections to the specific heat of the normal, degenerate component of 

■ quark matter. We study these non-perturbative non-Fermi liquid corrections 
(— i , to the specific heat and the temperature dependence of the chemical potential 
Qh! and show that these lead to a reduction of the specific heat. 

■ I. INTRODUCTION AND MOTIVATION 



A detailed description of the phase diagram of QCD in hot and dense environments 
can lead to a deeper understanding of Early Universe cosmology as well as the physics 
of compact stellar objects. Current experimental programs at CERN-SPS and BNL-RHIC 
probe QCD at high temperatures but relatively small baryon (and quark) density PJ and the 
future ALICE program at CERN-LHC will probe at even higher temperatures and smaller 
densities, in conditions that prevailed a few microseconds after the Big Bang. The goal of 
these ultrarelativistic heavy ion experiments is to search for the quark gluon plasma, a novel 
phase of QCD in which quarks and gluons are deconfmed. The region of the phase diagram 
of QCD for low temperatures T < 10 Mev but densities a few times larger than that of 
nuclear matter, cannot be studied with terrestrial accelerators but are the realm of compact 
stellar objects: neutron stars and pulsars. There is the tantalizing possibility that at the 
core of neutron stars where the density is up to 5 times that of nuclear matter, there could be 
a component of cold, degenerate quark matter with temperatures T < 1 Mev and chemical 
potentials y. ~ 300 - 500 Mev @,||. 
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An important aspect of dense QCD is the possibility of color superconductivity [|J- flT3"] . 
One gluon exchange leads to a pairing instability of the quark Fermi surface in the color 
antitriplet channel and while originally superconducting gaps were estimated to be of 
order A ~ 1 Mev more recent estimates including dynamical screening of the exchanged 
gluon || conclude that the gaps could be A ~ 50 — 100 Mev |5[]- [JHJ. 

The existence of quark matter phases in the cores of neutron stars and pulsars could 
have distinct observational consequences. The equation of state of quark matter may lead 
to pronounced delays in the spin-up history of neutron stars in low mass X-ray binaries 



and may explain recently observed anomalous frequency distributions ||14j| , a deconfinement 
transition can lead to observational consequences in the rotational properties of pulsars ||15| , 
superconducting quark matter can influence the strength and distribution of magnetic fields 
of pulsars anc ^ cou ld influence the cooling history of young neutron stars and 

protoneutron stars |19"1 , |20| . 



The properties of neutron stars and pulsars are studied by a satellite program that 
includes the Einstein Observatory, ROSAT, AXAF, RXTE and more recently CHANDRA 
and XMM that measure the (soft) X-ray emission from neutron stars, thus studying their 
cooling history as well as rotational properties (spin-up and spin-down) . 

A few seconds after the collapse of a Type II supernovae the newly born neutron star 
cools via the emission of neutrinos and antineutrinos and within few minutes the temperature 
falls to about 10 9 f^ (~ 0.1 Mev) cooling proceeds via neutrino emission for another 10 5-6 
years before photon emission becomes the most important cooling mechanism [ET ■ 



The most efficient neutrino emission mechanisms are the direct URCA processes fi + 1 — ► 
/2 + Vi ; f2^fi J rl J rVi with fi being either baryons or quarks and I either electrons or 



muons with v\ ; v\ their respective neutrinos and antineutrinos |LJ|. In the region between 
the crust and the core, where nuclear matter is the dominant component this process occurs 
with baryons and corresponds to the beta decay of the neutron and electron capture in a 
proton. If the core is composed of quark matter the corresponding processes are direct quark 
Urea, d — > u + e~ + v e ; u + e~ — > d + v e which have been studied in detail in ref. f2~8fl. The 



cooling rate is determined by the equation |2T|]- |26fl 



dU _ dT , r r . 

_ = c,- = -(U + 1,.) (l.l) 

with C v the specific heat (at constant volume and baryon number) and L vn are the neutrino 
and photon luminosities, any other potential sources of energy loss such as magnetic field 
decay, differential rotation etc. had been neglected in the above equation. 

The superconducting (and superfluid) component gives a contribution to the specific 
heat and the neutrino emissivity that is exponentially suppressed by the gap oc e _A//T if 
A ^> T [^,[M|,[n|. Recently, however, it was noted that novel aspects of cooling during the 
protoneutron star stage, when T ~ 30 — 50 Mev could be a result of the color superconducting 
transition [^] since in this regime A < T . 

For two light flavors (u,d) of quarks one gluon exchange leads to an instability in the 
antitriplet channel leading to a color superconducting state known as (2SC) [0,50]. This 
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state leaves quarks of one color ungapped while the others pair up with a gap A ~ 50 — 
100 Mev. The temperature of young neutron stars is T ~ 0.1 — lMev, hence in a color 
superconducting state the paired quarks will not contribute to cooling. Including the strange 
quark the color superconducting state for three flavors is a color-flavor-locked phase with all 
quarks paired. 

If the strange quark mass is too large to allow for S-wave pairing in the color-flavor-locked 



phase it may still form a P-wave condensate but with a very small gap |27j , in which case S- 
wave color superconductivity will be in the form of a (2SC) condensate. In this case the only 
color degree of freedom remaining "normal" and the strange quark contribute to the cooling 
via neutrino emissivity but the contribution from the strange quark will be suppressed by 
the (small) gap [27] as well as by the small Cabibbo angle [28[ . 



There are two important quantities that determine the cooling rate: the specific heat C v 
and the neutrino emissivity. Both are completely determined by the physics near the Fermi 
surface of the degenerate quarks and leptons (the only important ones are electrons). 



Iwamoto [fJHJ calculated the neutrino emissivity of normal quarks via the direct quark 
Urea processes, assuming the validity of Fermi liquid theory, and found it to be e oc T| 
(T 9 = T/10 9 K). The power of temperature is found on physical grounds: each degenerate 
fermion (quarks and electron) has associated a phase space factor of d 3 p/\p\ ~ ppT since only 
states within a region of width T near the Fermi surface contribute, this gives three powers 
of T, energy-momentum conservation restricts the neutrino to have energy-momentum ~ T 
and since the neutrinos are not Pauli blocked their available phase space is ~ E 2 dE u oc T 3 . 
The energy conservation delta function leads to a factor 1/T which is compensated by a 
factor E v ~ T in the numerator for the energy loss. Thus the energy loss by neutrino 
emission is oc Uf =Ui d, e {pF{f)T)T 3 ||28|| . The linear power of T arising from phase space for 
each degenerate fermion is determined by the single particle density of states near the Fermi 
surface which also determines the linear temperature dependence of the specific heat, as can 
be seen from the following argument. The contribution from massless degenerate quarks to 
the temperature change in the internal energy is determined by the states in a region of width 
« T from the Fermi surface, within which quarks have thermal energy « T. The density of 
single particle states at the Fermi surface pp = \i for relativistic fermions is oc pp = /i 2 , thus 
within a region of width T around the Fermi surface there are ~ /z 2 T states of energy ~ T 
leading to their contribution to the temperature change in internal energy 

5U q (/x, T) oc /i 2 T 2 -f C v oc n 2 T (1.2) 

The main reason for repeating these arguments is to highlight that both the neutrino 
emissivity and the specific heat of the normal component are determined by the (quasi) 
particle states within a region of width T near the Fermi surface. 

In ref. @J2^[ it was realized that the exchange of dynamically screened gluons that lead 
to the pairing instability also leads to a breakdown of the Fermi liquid picture for quarks 
near the Fermi surface. In particular these non-Fermi liquid corrections are responsible for 
a decrease in the gap fl2"9f . A detailed study of non- Fermi liquid aspects of the normal state 
was recently presented in P0[. There the spectral density, dispersion relation and width of 
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quasiparticles with momenta near the Fermi surface were derived at T = by implementing 
a renormalization group resummation of the leading logarithmic infrared divergences asso- 
ciated with the emission of soft dynamically screened transverse gluons. These non-Fermi 
liquid corrections modify the properties of quasiparticle states near the Fermi surface. Since 
cooling of degenerate quark matter is sensitive to the properties of the (quasi) particle states 
near the Fermi surface, it is natural to study the possible non- Fermi liquid corrections to the 
specific heat as well as the neutrino emissivity. 

Our goal in this article is to study in detail the non-Fermi liquid corrections to the 
specific heat to establish if these result in an enhancement or suppression of the cooling 
rate via neutrino emission through direct quark Urea processes, postponing the study of 
non-Fermi liquid effects upon the neutrino emissivity to a forthcoming article. 



II. QUASIPARTICLE DENSITY OF STATES, CHEMICAL POTENTIAL AND 

SPECIFIC HEAT: 

Our strategy to obtain the specific heat (at constant volume and baryon density) is to 
first obtain the internal energy and take its temperature derivative. The internal energy in 
turn is the expectation value of the full Hamiltonian H = H q + H g + H q _ g , the sum of the 
quark, gluon and interaction parts. The regime of interest for neutron star phenomena is 

T ~ 0.1 — 1 Mev ; fi q « 0.3 - 0.5 Gev (2.1) 



i.e, T <C fx. Furthermore the scale that enters in the resummed quark propagator [30] is 
M fa gn/27r at energy scales ~ 1 Gev the strong coupling constant a s ~ 0.6 hence M ~ fi q 
implying the hierarchy T <C M ~ fi. The contribution from massless gluons or Goldstone 
bosons remaining from the superconducting transition, to the internal energy is 

U g (T) oc T 4 (2.2) 

Therefore the gluon contribution to the internal energy is subleading by a factor T 2 /^ 2 
as compared to the quark contribution which, as argued above is 5U ~ fi 2 T 2 . Hence we 
will neglect the contribution from H g to the internal energy. The expectation value of 
H = H q + H q _ g can be obtained straightforwardly from the spectral density of the quark 
field by the following observation: the quark Hamiltonian including the interaction with the 
gluon field is bilineal in the quark field, and the equation of motion is 

therefore it follows that 

/ d 3 x^(x , t) 2 7 ° d t ty(x, t)=H (2.4) 

Hence we write 
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U(T, p) = (H) 



d 3 x idf 



x'=x\t=t' 



(2.5) 



The fermion correlator can be written in terms of the spectral representation of its spatial 
Fourier transform as follows l30|| 



-iS < {k;t-t') 
N(q ) 

Therefore, we find 
U(T,p) 



iSa,b( X X ) t 

dq p(q ,k)N(q )e- 1 ' 10 ^ 
1 



V 



U(T,p) 



d 3 k 



(2tt)* 



rfg goTr 7 p(q ,k) N(q ) 



In the grand canonical ensemble, the average number of particles is given by 



N = I d 3 x(^(f,t)7°^(£,£)) 



cfxTr 



-iS < {x-x , ;t-t')7° 



g=x';t=t' 

which is independent of temperature and given in terms of the spectral density by 

d 3 k 



N 



(27TY 



dq Tr p(q ,k)^ N(q 



(2.6) 
(2.7) 

(2.8) 
(2.9) 
(2.10) 
(2.11) 



The temperature dependence of the chemical potential p{T) is found by requiring that 
iV be independent of temperature. 

In the degenerate case temperature only affects states that are a distance ~ T from the 
Fermi surface, for positive baryochemical potential only quasiparticle states are important, 



since the contribution of the antiparticles will be suppressed by a factor exp — ^ <C 1 HlOlpT 
For massless quarks the spectral function p(qo, k) is of the form flQjll 

1 r 
2 



p(qo, k) 



V.(k)p.(q ,k) + V + (k) P+ (q ,k) 



(2.12) 



with V± = 7°±7-k and p± are the spectral density for quasi-antiparticles (+) and quasipar- 
ticles (-) respectively. As argued above temperature affects only states a distance ~ T from 
the Fermi surface, hence only quasiparticle states are important, and antiparticle states are 
exponentially suppressed. Therefore the contribution to the internal energy, specific heat, 
and (quasi) particle density for states near the Fermi surface are given by 

U(T, p) = 2j dq q ri(qo)N(q ) (2.13) 

C-% (2,4) 
N r 

A/"=- = 2 J dq oV (q )N(q ) (2.15) 



where we have introduced the single quasiparticle density of states 



d 3 k 



P-(?o,fc) 



(2.16) 



In perturbation theory, the quark self-energy features an infrared divergence for states 
near the Fermi energy which lead to the breakdown of the Fermi liquid picture. These 

divergences are a consequence of the exchange of dynamically screened soft transverse glu- 
ons, while longitudinal gluons are Debye screened. In reference |3(J the leading logarithmic 
infrared divergences were resummed via the renormalization group at T = 0, leading to a 
scaling form of the spectral density for frequency and momenta near the Fermi surface. The 
calculation for T 7^ but in the degenerate case T <C // is performed in a straightforward 
manner following the steps in this reference JJIJ. At finite temperature and chemical po- 



tential the self-energy for transverse gluons receives contributions from both the gluon and 
fermion thermal and dense loops [J31]J - PB"| . In the degenerate case derivatives of the Fermi- 
Dirac distribution function are sharply peaked at the Fermi surface. For N c colors and Nf 
quarks in the fundamental representation and qo, k ~ /j and after a renormalization group 
resummation of the leading logarithmic divergences along the same lines as that in |3(J we 
find that the improved quasiparticle spectral density for states near the Fermi surface is 
given by 



P-(Qo,k) 



sin[7rA] 






%\ 


90 

M 


-2A 


7T 


(% 


M 


-2A ~ \ 2 

— k cos[7rA] 1 -f 


(fcsin[7rA]) 



X-- 



■ Qo - // ; 

as N 2 C - 
6tt 2N c 
a. 



k = k 



M (T) = —Nf 
2ir 



/i 2 + — T 




(2.17) 

(2.18) 
(2.19) 

(2.20) 



where the Landau damping scale M 2 (T) displays the contribution from quark loops (first 
term) and gluon loops (second term) [0. This form of the spectral density is valid for 
l<7o|; 1^1 ^ M since the renormalization group resummation only includes the leading infrared 
logarithmic divergences and neglects perturbative analytic contributions. Since T <C M this 
form of the quasiparticle spectral density determines the finite temperature properties in the 
degenerate case. 

As discussed in detail in ref. 



30], this quasiparticle spectral density features a narrow 



resonance for k 7^ 0, at pO 



q (k) = sign(A;) \k\M~ zx cos[vrA] 



(2.21) 



Near the position of the resonance the spectral density can be approximated by a Breit- 



Wigner form 
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r iv 7 rri cos[7rA] 

P-(9o,K)l4b,fc«/x = Z p [k\ 



r(k) 



go(fc) 
M 



2 A 



Z ' [ * ] = (1 - 2A) 
r(jfc) = Z p {k}\k\sin[n\} 

Furthermore the group velocity near the Fermi surface is given by 

v g (k) = -Jr = Z p [k] cos[ttA] 



(2.22) 

(2.23) 
(2.24) 

(2.25) 



Therefore the residue of the "quasiparticle pole", the "quasiparticle width" and the 
"quasiparticle group velocity" vanish near the Fermi surface as [ |30[ 



2A 
-2A 



Z p [k] oc \k — /cf| t 
r(fc) oc | A; — k F \ T =^ 
v g (k) = Z p [k] cos[7rA] oc \k — /cp| T 



2A 

-2A 



(2.26) 
(2.27) 
(2.28) 



The vanishing of the residue at the quasiparticle pole at the Fermi surface is the hallmark 
of the breakdown of Fermi liquid theory. In a degenerate Fermi gas, all of the thermodynamic 
response functions are determined by the quasiparticle properties near the Fermi surface. The 
vanishing of the quasiparticle residue at the Fermi surface, and the consequent breakdown of 
Fermi liquid theory has the potential for modifying the thermodynamic response functions. 
The important aspect that we seek to understand is how these non-Fermi liquid features 
affect the specific heat of quark matter which is relevant for the cooling of neutron stars 
with a normal component of quark matter in the core. 

As discussed in |$0| the non-Fermi liquid corrections to the quasiparticle spectral den- 
sity are non-perturbative as explicitly displayed by the anomalous powers in the dispersion 
relation and the wave function renormalization. The renormalization group improved quasi- 



particle spectral density ( 2.17 ) has been obtained by a resummation of the leading infrared 
divergences near the Fermi surface and neglects perturbative corrections which are analytic in 
the frequency and momentum near the Fermi surface |30J . The validity of ( 2.17|) is restricted 
to a region — M < qo', k < M where the contribution from the dressed transverse gluons 
is dominated by Landau damping 3JJ. Away from the Fermi surface the spectral density 
must match smoothly to that of the fermionic collective modes (see [j30l , |33 -35l) which for 
frequency and momenta » M lead to perturbative corrections to the quasiparticle spectral 
density for states with frequency and momenta away from the Fermi surface. 

Thus we will restrict our study of the specific heat to the non-perturbative aspects of the 
non-Fermi liquid corrections by focusing on the contribution to the specific heat of states 
within a width of order M near the Fermi surface and on the terms with anomalous ex- 
ponents which are the hallmark of non-Fermi liquid corrections. Therefore, we will neglect 
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the contributions from perturbative terms since these will arise from: i) perturbative cor- 
rections to the quasiparticle spectral density not included in ([2.17D and ii) the contribution 
from states a distance > M below the Fermi surface, none of which are associated with the 
breakdown of the Fermi liquid description. 



A. Single quasiparticle density of states 



The expressions for the internal energy, specific heat and density given by ( p. 13 - 2. 15 ) 
require the computation of the single quasiparticle density of states (]2.16|) . The non- Fermi 
liquid features of the single quasiparticle density of states can be extracted by restricting the 
momentum integral to states within a width of order M around the Fermi surface, hence 
non-Fermi liquid contributions to i](qo) are given by 



v(Qo) 



2tt 2 



ai 



dk 



M 



~l 2 



H + h p-(q ,k) 



(2.29) 



Straightforward integration leads to the following explicit form 



MY 

71 J 
1 



+ 



2 shi7rA 



sin7rA f f u . 
{ z — z \ —a + z cos 7tA ArgTh 

7T I \M ' 



M) 



2z— acosrcX + z cos27rA 
M 



2 Z COS7tA" 

1 + z 2 
arctan 



2 z sin 7rA 



7T 



(2.30) 



with 



-j^-j ; a = sign(g ). 



1-2A 



Since the specific heat is determined by a small region of width q 
imate the above form of f](q~o) for z 1 



(2.31) 

T<cMwe approx- 



7T 

M 

^2 



M 



7T Z 



M 



1 — sin 7rA + \i a cos tt\ 



M 2 



7T 



2ua M 

sm 27rA cos 27rA 

7T 2 



0{z 2 



(2.32) 



This expression has the correct free-field limit, which corresponds to p~{qo, k) = 5(qo — k). 
A noteworthy point to emphasize is that while the quasiparticle residue vanishes at the 
Fermi surface, the single quasiparticle density of states is non-vanishing and finite but with 
a divergent frequency derivative at the Fermi surface (go = 0). This divergence in the 
frequency derivative introduces non analytic finite temperature corrections to the chemical 
potential which are explored below. 
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B. Chemical Potential: 



In the grand canonical ensemble, the temperature dependence of the chemical potential 
is obtained by requiring that the mean number of particles is fixed and is independent of 
temperature. The chemical potential is the Lagrange multiplier that enforces the constraint 
of a given particle number and its temperature dependence is obtained from the condition 



dM 
~dT 



= J dq 



dN(qo) drj(q ) 



(2.33) 



with Af is the particle density given by (|2.15| ). 

In the degenerate case dN(q )/dT is localized in a region of width ~ T around q = \i 
hence the first term in ( |2.33| ) above will only receive contributions from quasiparticles near 
the Fermi surface for which the single quasiparticle density of states is given by (|2.30| ). 

While in the free theory the single particle density of states does not depend on the 
temperature or chemical potential, in the interacting theory it depends on T both through 
M(T) given by ( |2.2(J| ) as well as through the chemical potential, thus the second term in 
( |2.33| ) is non- vanishing. However, while the second term in ( |2.33|) receives contributions from 
states well below the Fermi surface, the non-perturbative non-Fermi liquid corrections are 
only associated with states within a region of width M near the Fermi surface. Quasiparticle 
states or collective modes well below the Fermi surface give perturbative corrections which 
are analytic in q Q in the region %/M <C 1. Such contributions will be neglected in this 
article. 

The contribution to drj/dT from the quasiparticle states near the Fermi surface given by 
eq. (|2.30| ) can be obtained as follows. Writing 



dr](q ) _ dri(q ) dp drj(q ) dM 
dT dp dT dM dT 1 ' ' 

it is readily apparent from eq. ( |2.30| ) or its small z approximation eq.( |2.32| ) that 

dy(qo) . dM T 

— r— — oc A ; — — oc A— 2.35 
dM ' dT M v ' 

therefore the second contribution in eq. (|2.34j) is perturbatively small 0(X 2 T/M), does not 
lead to anomalous powers of the ratio T/M < 1 and will be neglected, consistently with 
neglecting perturbative corrections in the spectral density and the quasiparticle density of 
states. Upon integrating in q the term N(q~o) up to a cutoff of order ~ M, one finds 
that the contribution from the first term of ( 2.34Q is again perturbative, and of the same order 



as terms that have been neglected in the spectral and quasiparticle densities. Therefore, we 
conclude that the contribution of the second term in eq. ( |2.33| ) is of the same order as terms 
that were neglected in the perturbative expansion and hence it will be consistently neglected. 

Hence, we are led to the following temperature derivative of the chemical potential to 
leading order 
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dn(T) _ J- M dq Q r](q )D(q )(Pq ) 



dT 



; D{x) 



(2.36) 



The leading contribution in the weak coupling and degenerate limit can be extracted by 
changing variables (3q~o = x and taking M/T — > oo. We find 



±, ( T) = -4 C08 ,A F(A) (-^) ( £ + O U-lJ (£ 



-4A X 



(2.37) 



leading to the following form of the temperature dependence of the chemical potential 



AT) = Ao) 



T 



—2 A 



+ o 



T 



Xo)J (m(o)J 



-4 AN 



where we have introduced the functions 

F(X) = T(3 - 2 A) C(2 - 2 A) (1 - 2 2A " X ) 

k(A) = 8 cos(ttA) T(2 - 2 A) C(2 - 2 A) (1 - 2 2X - 1 ] 



4 COS7TA 

1 - A 



F{\) 



(2.38) 

(2.39) 
(2.40) 



in terms of the Gamma and Riemann zeta functions. 

The non-analiticity of the temperature derivative of the chemical potential reflects the 
divergence of the frequency derivative of the single quasiparticle density of states near the 
Fermi surface. 



C. Specific Heat: 



The specific heat (per unit volume) given by ( |2.14j ) can be written as 
C v = 2 J , 



^ dN(q ) d7](q ) 
r]{q )——- + — N(q ) 



dT 



dT 



(2.41) 



To extract the leading order in non-Fermi liquid corrections we follow the analysis invoked 
in obtaining the temperature derivative of the chemical potential. In the degenerate case the 
first term in eq. (|2.41|) above is determined by the quasiparticle density of states eq. ( p.30|) 
for % ~ T M. An estimate of the second term is obtained by following the arguments 
presented for the case of the chemical potential. A similar analysis reveals that the second 
term is perturbative in the weak coupling limit and of the same order as terms that were 
neglected in obtaining the spectral density and the single quasiparticle density of states near 
the Fermi surface. Therefore the second term in eq. ( |2.41| ) will be consistently neglected. 

Using the condition ( |2.33j ) and the fact that near the Fermi surface the quasiparticle 
density of states is a function of %, the leading order (in coupling) contribution to the 
specific heat is obtained from 
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a 



M 

- - , dq q ri(q )f3D(q ) 

M 



Pqo + 



dfi(T) 
dT 



(2.42) 



with D(x) given in eqn. ( |2.36| ). 

After a rescaling q = xT, using the result (|2.37| ) for d ^jp and neglecting terms pertur- 
bative in the coupling A we find 



1 + 



6 

V 2 

2 



F(2 A - 1) cos 2ttA + 8 F(A) 2 cos 2 ttA 



' t y ( _r_ \ 



-4A 



(2.43) 



with the functions F(X) ; k(X) given by eq. ( |2.39|J2.40| ) above. 

The non-Fermi liquid contributions are explicit in the anomalous power laws with the 
temperature, which are a result of the divergent frequency derivative of the single quasipar- 
ticle density of states near the Fermi surface. We emphasize that for T<M the non-fermi 
liquid corrections to the specific heat are much larger than those of the massless gluons or 
any other Goldstone bosons to the specific heat, which is of order T 3 . 

To illustrate the magnitude of the non-Fermi liquid corrections to the specific heat we 
show in fig. 1 the function G(z) = 3C V / fi 2 (0)T vs. z = T//i(0) for three colors and two 
flavors of massless quarks. We have used the one loop QCD running coupling at a scale 
/i(0) without further justification^]. The solid line corresponds to /x(0) = 0.5 Gev and the 
dashed line to /i(0) = 1 Gev, for these values of the strong coupling constant, the value 
of the effective coupling at these scales are A = 0.152 for \x = 0.5 Gev and A = 0.046 
for fi = 1 Gev therefore the reliability of the perturbative calculation for the values of the 
chemical potentials expected at the core of neutron stars is at best questionable. This issue 
notwithstanding this qualitative estimate points out that non-Fermi liquid corrections to 
the specific heat result in a change of a few percent. The physical reason for this result 
is that whereas the quasiparticle residue (wave function renormalization) vanishes at the 
Fermi surface, the single quasiparticle density of states is non-vanishing and finite, but its 
frequency derivative at the Fermi surface is divergent. The finite, non-zero value of the 
density of states implies that the non-Fermi liquid features will only result in corrections to 
the Fermi liquid specific heat. The numerical estimate, with the caveat that its regime of 
validity is restricted to small coupling, indicates that these corrections are rather small. 

Keeping higher order terms in the perturbative expansion, will not dramatically modify 
this conclusion since our analysis captures the most relevant aspects of the non- Fermi liquid 
corrections. For the values of temperature and chemical potential corresponding to young 



1 Reference [37] gives a thorough discussion of the dependence of the coupling on the chemical 
potential and the inherent ambiguities in the choice of scale. 
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neutron stars we see that the non-fermi liquid corrections amount to a few percent, unlikely to 
make a dramatic change in the cooling history of young neutron stars. For protoneutron stars 
with higher temperatures T ~ 30 — 50 Mev and smaller chemical potentials, the corrections 
could be of order 10 — 15% but in this regime the reliability of the calculation is questionable. 
We emphasize that this numerical analysis is only intended to provide a qualitative estimate 
of the magnitude of the non- Fermi liquid corrections but should not be taken as a quantitative 
indicator since higher order terms, which are subleading and have been neglected in the 
perturbative expansion can become comparable in this range of couplings. 



D. Discussion 



The fact that the wave function renormalization, i.e, the residue at the "quasiparticle 
pole" vanishes at the Fermi surface as given by eqn. Q2.26Q would seem to indicate that 
quasiparticle density of states should vanish near the Fermi surface with an anomalous 
power law. This in turn would suggest that the linear power of temperature, which is the 
leading contribution in the degenerate limit, would be multiplied by an anomalous power 
of T I /i. However, the reason that the quasiparticle density of states remains finite at the 
Fermi surface despite the vanishing of the quasiparticle pole is because the group velocity 
vanishes proportionally to the wave function renormalization as displayed by eqn. ( |2.28| ). 
To see this feature clearly, we can approximate the spectral density ( [2.22|) in the limit of 
vanishing width by 



fc^O 



Z p [k] cos[7i\]5(q 



&(*)) 



(2.44) 



with %{k) given by eqn. ( |2.2ip . The integral over k in eqn. ( |2.29| ) can be performed by 
changing variables to %{k) where the Jacobian of the transformation is precisely l/v g with 
v g the group velocity and by eqn. ( 2.2S ) this Jacobian exactly cancels the wave-function 



renormalization. This analysis clearly reveals that the reason that the quasiparticle density 
of states at the Fermi surface remains finite is the cancellation between the residue at the 
quasiparticle pole and the group velocity, which both vanish proportionally near the Fermi 
surface. 

At this stage it is also important to compare our results to those obtained in reference |RS 



in which the infrared divergences associated with magnetic photon exchange were studied 
for the first time within the context of the non-relativistic electron gas. There are several 
important differences between our treatment and results and those of reference ||38|| . i) 
In reference [[38|] the temperature derivative of the chemical potential was neglected in the 
derivative of the grand potential with respect to temperature (see equations (5.6) and (5.7) 
of reference |38|]). This is an important source of discrepancy because the temperature 
derivative of the chemical potential has a linear term in the temperature that cancels against 
a similar one arising from derivatives of the density of states and the distribution function. 
We have included the temperature derivative of the chemical potential in the computation of 
the specific heat explicitly through the condition ( |2.33|) . ii) The calculation in |38] missed the 
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fact that the vanishing of the group velocity is compensated by the vanishing of the residue 
at the quasiparticle pole as explicitly discussed above, iii) While in reference |38j only the 
perturbative logarithmic terms were accounted for, the renormalization group resummation 
makes clear that these logarithms sum up to anomalous power laws, a naive expansion in 
terms of the coupling constant reveals those logarithms, but the resummation makes explicit 
the vanishing of both the residue and the group velocity. 

A consistent resummation and analysis of the quasiparticle pole and group velocity as 
explicitly discussed here reveals that the density of states remains finite at the Fermi surface 
and the non-Fermi liquid corrections are subleading. 



III. CONCLUSIONS 

The exchange of dynamically screened gluons leads to a pairing instability and at the 
same time to a breakdown of the Fermi liquid description of the normal state of QCD. 
The superconducting component does not contribute to the cooling of young neutron stars 
because of the exponential suppression of the specific heat and the neutrino emissivity. If 
color superconductivity is in the form of a 2SC condensate, or some quarks do not pair, 
the normal quarks contribute to cooling of the neutron star via neutrino and antineutrino 
emission. 

The search for novel observational aspects of quark matter at the core of neutron stars 
motivated us to study the effects of the breakdown of the Fermi liquid description of quasi- 
particle states of the normal component near the Fermi surface in the specific heat. This 
aspect is relevant to understand the cooling rate through direct quark Urea processes. Our 
analysis is based on a renormalization group improvement of the single quasiparticle spectral 
density in the degenerate case and reveals that non-Fermi liquid corrections tend to decrease 
the specific heat and therefore increase the cooling rate. However these corrections are rather 
small, about 10 — 15% even in the most favorable scenarios, and probably unobservable in 
the cooling history of neutron stars or even protoneutron stars. A complete assessment of 
potential non-Fermi liquid corrections to the cooling rate still requires an understanding of 
these corrections on the neutrino emissivity. Both quarks and leptons (electrons) that enter 
in the quark direct Urea process will receive non- Fermi liquid corrections through their spec- 
tral density although the corrections for quarks, in terms of the strong coupling constant are 
more important than those for the electrons. We are currently studying these corrections in 
the weak matrix elements for the quark direct Urea processes and expect to report on our 
results in the near future. 
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